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ABSTRACT 

We present the first results of a search for Lya emitters (LAEs) in the DEEP2 spectroscopic database 
that uses a search technique that is different from but complementary to traditional narrowband 
imaging surveys. We have visually inspected ~20% of the available DEEP2 spectroscopic data and 
have found nine high-quality LAEs with clearly asymmetric line profiles and an additional ten objects 
of lower quality, some of which may also be LAEs. Our survey is most sensitive to LAEs at z=4.4-4.9 
and that is indeed where all but one of our high-quality objects are found. We find the number 
density of our spectroscopically-discovered LAEs to be consistent with those found in narrowband 
imaging searches. The combined, averaged spectrum of our nine high-quality objects is well fit by a 
two-component model, with a second, lower-amplitude component redshifted by ^420 kms -1 with 
respect to the primary Lya line, consistent with large-scale outflows from these objects. We conclude 
by discussing the advantages and future prospects of blank-sky spectroscopic surveys for high- z LAEs. 
Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — techniques: 
spectroscopic 



1. INTRODUCTION 

Over the last decade, studies of how galaxies form 
and evolve have been pushed to unprecedented look-back 
times, allowing us to observe galaxies at epochs when the 
universe was only a fraction of its present age. A number 
of surveys for high-redshift (z>2) galaxies has been un- 
dertaken, including many in the rest-frame UV/optical 
(e.g., Steidel et al. 1996, 1999, 2003; Sawicki et al. 1997; 
Franx et al. 2003; Daddi et al. 2004; Sawicki & Thomp- 
son 2005; Iwata et al. 2007) and far-IR (e.g., Barger et al. 
1998; Hughes et al. 1998; Blain et al. 1999; Smail et al. 
2002; Chapman, et al. 2003; Webb et al. 2003; Sawicki & 
Webb, 2005; Coppin et al. 2006). The study of galaxies 
selected by their strong Lya emission allows an investiga- 
tion of a population that's potentially very different from 
these other high-z galaxy populations. Lya selection not 
only gives us the prospect of probing galaxies with very 
faint spectral continuum levels, but also may select sys- 
tems in very early stages of a starburst (e.g., Malhotra & 
Rhoads 2002). It is thus not surprising that over the last 
several years much observational effort has gone into Lya 
surveys (e.g., Ajiki et al. 2003; Hu et al. 1998, 2004; Ku- 
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dritzki et al. 2000; Ouchi et al. 2003; Rhoads et al. 2000, 
2003; Rhoads & Malhotra 2001; Santos et al. 2004). In- 
deed, the short rest-frame wavelength of the Lyal216 
line allows galaxies to be detected and spectroscopically 
confirmed to very high redshifts, z<7, without the need 
for infrared observations. 

The ability to push to such high redshifts raises the 
intriguing possibility of probing the tail end of the epoch 
of reionization, regarded by many as marking the onset 
of galaxy formation in the Universe. The expected evolu- 
tion of the IGM neutral fraction across the reionization 
epoch has led to the prediction that the observed Lya 
luminosity function (LF) will undergo rapid evolution as 
IGM opacity to Lya photons drops with decreasing red- 
shift (e.g., Haiman & Spaans 1999; Malhotra & Rhoads 
2004; Kashihawa et al. 2006) and attempts have been 
made to detect the signature of this effect albeit with con- 
tradictory results (Malhotra & Rhoads 2004; Kashikawa 
2006). At slightly lower redshifts, the LF of high- z galax- 
ies may be evolving due to evolution in the intrinsic prop- 
erties of galaxies (e.g., Sawicki & Thompson 2006; Iwata 
et al. 2007; see also Bouwens et al. 2007). Consequently 
even if the observed evolution of the Lya LF proves to 
be real, it will remain an open question as to whether it 
reflects an evolution in the opacity of the IGM, or evo- 
lution in the intrinsic properties of galaxies at a time 
when they were young. To believe the former we must 
understand the latter, and thus it is important to study 
the properties of Lyman Alpha Emitters (LAEs) in de- 
tail not only at very high redshift, z>6, but also as a 
function of time. 

To date, narrowband imaging searches for high-z LAEs 
have been the most fruitful approach to finding LAEs 
and have yielded significant samples of galaxies at z~3- 
6.5, with several dozen objects now confirmed spectro- 
scopically by various teams (e.g., Rhoads et al. 2003; Hu 
et al. 2004; Dawson et al. 2004, 2007; Shimasaku et al. 
2006; Kashikawa et al. 2006). An alternative to narrow- 
band imaging surveys lies in direct, blank-sky spectro- 
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scopic searches whose strategy is to bypass the photome- 
try stage and find galaxies directly in blank-sky spectra. 
One advantage of this approach is the potential saving of 
the time associated with deep narrowband imaging; an- 
other is the gain in sensitivity that results from working 
at the natural linewidth of the line (~10-20A) rather 
than fighting against the sky background admitted by 
the full width of the narrowband filter (~100-200 A) (see 
Martin & Sawicki 2004). 

Until recently, the chief disadvantage of spectroscopic 
searches has been the small sky area accessible to most 
spectrographs. Serendipitous spectroscopic discoveries 
of single emission-line objects (e.g., Franx ct al. 1997; 
Dawson et al. 2002), while useful, do not yield sam- 
ples of galaxies suitable for statistical analysis. Spec- 
troscopic searches of regions enhanced by gravitational 
lcnsing have met with some success (Santos et al. 2004; 
Stark et al. 2007), probe very deep into the luminosity 
function but over very small areas. Until recently, mul- 
tislit blank-sky searches (e.g., Crampton & Lilly 1999; 
Martin & Sawicki 2004; Tran et al. 2004) have found 
many line emitters but no confirmed LAEs; only now are 
the sky areas accessible to this technique becoming large 
enough to report successful detections (Martin, Sawicki, 
Dressier, & McCarthy 2008). 

In the present work we take an approach that is hybrid 
to the dedicated blank-sky spectroscopic surveys and the 
serendipitous discoveries. Specifically, we carry out a 
dedicated search for serendipitous LAEs in the extremely 
large, existing spectroscopic database of the DEEP2 sur- 
vey. The present paper reports our first discovery of a 
sample of Lya emitters using these data and illustrates 
some of the constraints that this approach can place on 
the properties of these objects. The paper is organized 
as follows. In §2 we describe our DEEP2 data and their 
suitability to the task at hand. In §3 we describe our 
search for LAEs and report on the objects that we have 
found. In §4 we present a number of tests with which 
we build confidence in the high-redshift nature of our 
objects. In §5 we constrain the number density and lu- 
minosity function of high-z LAEs and in §6 we study 
the spectral properties of the Lya emission line at high 
redshift. In §7 we discuss future prospects for finding 
LAEs using deep blank-field spectroscopy and summa- 
rize our main findings. Throughout this paper we adopt 
the 12m=0.3, flA—0.7, H o —70 km s" 1 Mpc" 1 cosmology. 

2. DATA 

Our study is based on an existing and very large 
database of two-dimensional spectra that were collected 
as part of the DEEP2 redshift survey. We searched 
these spectra to identify emission-line objects that are 
not (other than by proximity in projection) associated 
with these primary DEEP2 spectroscopic targets. In this 
section we first describe the DEEP2 survey itself before 
launching into description of our search for serendipitous 
high-redshift galaxies in §3. 

2.1. The DEEP2 Survey 

The DEEP2 survey (Davis et al. 2003, 2005; Coil et al. 
2004) is a large spectroscopic redshift survey designed to 
measure redshifts and other properties of ~50,000 galax- 
ies using DEIMOS (DEcp Imaging Multi-Object Spectro- 
graph; Faber et al. 2003) on the 10m Keck II telescope. 



The primary goal of DEEP2 is the study of galaxy evo- 
lution at intermediate redshifts, z~l. 

The survey spans four fields, with Field 1 (the Ex- 
tended Groth Strip) being a strip of 0.5 deg 2 and Fields 
2, 3, and 4 each spanning 1 deg 2 , although in Fields 2 
and 4 spectroscopic coverage was completed over only 
^0.67 deg 2 . Primary targets for spectroscopy are se- 
lected in deep BRI images taken with the CFH12K cam- 
era on the Canada- France-Hawaii Telescope (for details, 
see Coil et al. 2004, Davis et al. 2005). These images are 
used to select galaxies for DEIMOS spectroscopy down to 
a limiting magnitude of Rab — 24.1. Except in Field 1, 
a two-color cut is also applied to exclude galaxies with 
redshifts z<0.75. 

Spectroscopy is obtained in multislit mode through 
custom-designed slitmasks. Slits arc 1.0" wide and 
around 140 slits are cut per DEIMOS mask. Slit lengths 
range from 2 to 67 arcsec, with 90% of slits being between 
3.7 and 11.3 arcsec long and only 0.4% of them shorter 
than 3 arcseconds. Typical exposure times are 1 hour, 
although, when observing conditions are sub-optimal, 
this exposure time is usually increased to achieve a con- 
sistent survey depth across all masks. The primary 
intermediate-redshift target galaxy in each slit occupies 
only part of the slit length, and substantial area in each 
slit is exposed to blank sky in order to allow for accurate 
sky subtraction in the spectral reduction process. In our 
work this "background" blank sky becomes the primary 
resource as it is there that we search for high-redshift 
emission-line galaxies. 

The DEEP2 survey uses a 1200 line mm" 1 grating, 
which covers a spectral range of 6400-9000 A at the mod- 
erately high resolution of R=5000. The spectroscopic 
data were pipeline-reduced using an IDL pipeline de- 
veloped at UC-Berkeley (Cooper et al. 2008, in prepa- 
ration 9 ). For our emission-line search we use the data 
that have been taken through the reduction process that 
includes standard CCD pre-processing followed by two- 
dimensional sky subtraction. The sky subtraction pro- 
cess (S. Buries and D. Schlegel 2008, in preparation) 
works exceptionally well and brings the data very close to 
the Poisson sky limit, with very few sky-subtraction arti- 
facts remaining. It is these 2-dimensional sky-subtracted 
spectra that we use to search for high-z line-emitting 
galaxies. 

2.2. DEEP2 as a high-z Lya survey 

While the chief goal of the DEEP2 survey is the study 
of the assembly of galaxies at intermediate redshifts, its 
extensive wavelength coverage, high spectral resolution, 
and large area make it an excellent dataset to search for 
galaxies at much higher redshifts. 

The DEEP2 spectral range of -6400-9000 A allows 
us to search for LAEs in the redshift range z~4.2-6.6. 
The large wavelength coverage of each spectrum also al- 
lows us to immediately reject most foreground interlop- 
ers through the detection of associated spectral lines at 
other wavelengths. Moreover, the spectral resolution of 
the data is high enough to resolve line shapes and, in 
fact, in many cases our candidates exhibit the asymmet- 

9 For a description of the DEEP2 spectroscopic pipeline see 
http : //astro . berkeley . edu/~cooper/deep/spec2d/primer . html 
and http : //astro . berkeley . edu/~cooper/deep/spec ld/primer . html 
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ric shape expected of the Lya line. Indeed, at our res- 
olution of R=5000 we can easily resolve the [O II] 3727 
doublet that can be a major source of contamination in 
Lya surveys (see, e.g., Martin & Sawicki 2004; Martin, 
Sawicki, Dressier, & McCarthy 2006, 2008); while the 
doublet may sometimes be smeared due to galaxy inter- 
nal kinematics, our high spectral resolution allows us to 
reject many, if not all, of the contaminating [O II]3727 
emitters. Finally, the availability of BRI imaging allows 
us a further check on our candidates since any bona-fide 
z>4 emitters must also be B-band drop-outs. In contrast 
to most LAE surveys, the DEEP2 spectroscopic data are 
already in hand and allow us to perform all these rejec- 
tion tests without the need for additional observations. 

The DEEP2 survey has accumulated a very large slit 
area exposed to the sky, which, combined with its wide 
wavelength coverage, results in a very large LAE survey 
volume. Each DEEP2 DEIMOS mask contains around 
140 l"-wide slitlets, an arrangement that, after account- 
ing for DEIMOS' inter-chip gaps — gives the same sky 
coverage as a 15.4 arcminute-long l"-wide longslit - 
namely 0.257 arcmin 2 per mask. The full DEEP2 sur- 
vey consists of 385 masks, giving a sky coverage of 98.8 
arcmin 2 . The present study is based on 83 masks, which 
corresponds to 21.6% of that 385-mask total. Most of 
the spectra surveyed in the present work come from 
Field 4 of the DEEP2 survey. The area surveyed in the 
present work is 21.6 arcmin 2 , which is somewhat larger 
than the areas surveyed by the recent custom-designed 
LAE multislit spectroscopic searches of Martin & Sawicki 
(2004; 5.1 arcmin 2 ) and Tran ct al. (2004; 17.6 arcmin 2 ), 
but substantially smaller than the —200 arcmin 2 of the 
current state-of-the-art multislit search of Martin et al. 
(2008). However, while the dedicated multislit surveys 
typically target only a narrow spectral range of 100- 
200 A, our search is sensitive to objects over —6400- 
9000 A. Consequently, our search spans a disproportion- 
ately larger volume than is suggested by a simple com- 
parison of survey areas. Our search — which uses only 
—20% of the available DEEP2 data — surveys a volume 
of 6.9 xlO 4 Mpc 3 and is thus already comparable to the 
4.5xl0 4 Mpc 3 of the state-of-the-art dedicated Magel- 
lan multislit survey of Martin et al. (2008). It is much 
bigger than the volumes covered by the earlier dedicated 
multislit searches. 

The top panel of Figure 1 shows the volume over which 
the DEEP2 survey is sensitive to LAEs. This volume 
was calculated from a representative sampling of actual 
DEEP2 masks by multiplying the spectral coverage of 
each slit in these masks by that slit's length and width, 
and converting the result into volume as function of wave- 
length. The upper curve shows the volume for the entire 
(385-mask) DEEP2 survey and the lower curve that of 
the present search that uses only 83 of the DEEP2 masks. 
The gray vertical lines show the locations of bright night- 
sky emission lines. 

Several interesting features can be seen in Fig. 1(a), 
and here we briefly explain their origins. The steeply 
sloping ends of the wavelength range are a result of sum- 
ming over the many thousands of spectra that make up 
the survey: while each individual DEEP2 spectrum has a 
well-defined minimum and maximum wavelength, these 
wavelengths vary from spectrum to spectrum subject to 
the slit's location on the mask and hence the spectrum's 



position with respect to the edges of the DEIMOS CCDs. 
Taken as an average over the survey, the resulting sur- 
vey volume has endpoints that are smoothly declining 
rather than abrupt as would be the case for a single spec- 
trum. Next, the smooth, shallow decline of the survey 
volume per redshift interval over the bulk of the survey 
range (z=4.6-6) is primarily a reflection of the redshift 
dependence of comoving volume per unit sky area in the 
adopted cosmology. And finally, the shallow dip in cov- 
erage at —7800 A is a consequence of the inter-chip gap 
between the two rows of DEIMOS CCDs. As with the 
spectral endpoints discussed earlier, this coverage gap 
is a sharply defined wavelength region for each DEEP2 
spectrum, but taken as a sum over all the DEEP2 spectra 
this feature becomes a broad, shallow dip. 

The vertical gray lines in Fig. 1(a) mark regions of in- 
creased noise due to night sky lines. In most surveys, 
these sky lines would make most of the wavelength range 
unusable for line emitter searches and this is why most 
LAE surveys target the few spectral regions that are rel- 
atively skyline- free (e.g., —8200 A, —9200 A). However, 
the DEEP2 spectral resolution is sufficiently high that in 
our case the skylines contaminate only a fraction of the 
total wavelength range and we are able to search for line 
emitters in the many inter-skyline regions throughout 
—6400-9000 A. Figure 1(b) shows the cumulative volume 
of the present, 83-mask survey (integrated from lower to- 
wards higher redshift), with the dashed line representing 
the volume as it would be in the absence of skylines and 
the solid curve showing the volume once we assume that 
the regions of high sky noise are rendered unusable by 
the higher sky noise. At the DEEP2 resolution, sub- 
stantial useful volume remains even in the sky-line rich 
regions. Having said that, however, we note that the 
region <7200A is particularly skyline-free, resulting in 
a concentration of useful survey volume at z— 4.6. It is 
at these redshifts that most of our highest-quality LAE 
candidates are found. 

Overall, the DEEP2 survey provides a very impressive 
and unique database for line emitter searches. It covers 
a large volume, giving us the prospect of finding sta- 
tistically significant samples of objects. It has a large 
wavelength baseline and thereby offers the potential for 
evolutionary studies. It has deep BRI imaging that helps 
in the rejection of foreground interlopers. It reaches faint 
spectroscopic flux limits of a few xlO -18 ergcm _2 s _1 , 
as we illustrate in § 3. And the moderately high spec- 
tral resolution of its data not only provides us with large 
skyline-free regions to search in, but also gives us the 
ability to immediately ascertain the identity of our LAE 
candidates by inspecting their line profiles without the 
need for additional, follow-up spectroscopy. 

3. LYMAN ALPHA EMITTERS IN DEEP2 

3.1. The search 

We searched the two-dimensional sky-subtracted 
DEEP2 spectra for serendipitous line-emitting objects. 
The first steps of the search, performed by one of us (BL) 
were to visually identify an initial sample of serendipitous 
emission lines and then to define a preliminary catalog 
of those most likely to be LAEs on the basis of a number 
of interloper-rejection criteria. Subsequently, this pre- 
liminary catalog was carefully re-examined by three of 
us (BL, MS, PG) to arrive at a final, robustly conserva- 
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tive LAE catalog. Our guiding principle has been to re- 
ject objects liberally rather than to retain objects about 
whose nature there was some question. This philosophy 
results in an underestimate of the number of true LAEs 
but this is a deliberate policy on our part. The steps of 
our classification procedure are explained in the following 
paragraphs. 

The first and most time-consuming step (~400 person- 
hours) in the search was the visual inspection of over 
10,000 two-dimensional spectra, representing ^20% of 
the DEEP2 dataset. This search focused on identify- 
ing emission lines that were clearly not associated with 
the primary spectroscopic DEEP2 targets. The emission 
lines of interest here were either lines that were offset 
spatially along the slit from the primary target, or lines 
that were superimposed on the spectrum of the primary 
target but clearly did not match any lines consistent with 
the redshift of the primary. 

Once a serendipitous candidate had been identified as 
described above, a number of tests were applied to ascer- 
tain its identity. First, we rejected from our LAE catalog 
any candidate (other than those few that were spatially 
coincident with foreground objects) that had one or more 
additional emission lines at the same spatial location in 
the 2D spectrum: for a true LAE galaxy the DEEP2 
wavelength coverage would not include any prominent 
lines other than the Lya, and thus any lines additional to 
the discovery line brand the candidate as a low- z object. 
The rejection procedure involved both the inspection of 
the entire length of the 2D spectrum and a targeted in- 
spection of special regions in the 2D and ID spectra. 
These special regions were regions where other emission 
lines would be expected to lie if our candidate were not 
Lya, but, rather, [O III] at 5007 A, H/3 at 4861 A, Ha at 
6563 A, or the [O II] doublet at 3727 A. We also verified 
that the candidate is not an [O II] emitter by checking 
for the double-peaked profile associated with this doublet 
line. In some cases, motions within galaxies can smear 
the doublet, but for all but very high electron densities, 
smeared galactic [O II] can be expected to have a profile 
that is asymmetric in the opposite sense to that of Lya 
because of the relative intensities of the [O II]3726,3729 
doublet components (e.g., Osterbrock 1989); such in- 
verted [O II] ratios are found to be very rare observa- 
tionally (Weiner et al. 2006). Thus, to check for smeared 
[O II], we verified that our candidates do not have the 
[O II]-like blue-tailed line profiles and found none. Next, 
we also checked that neither the 2D nor ID spectra had 
continuum blueward of the candidate line. Once the can- 
didate met all of these spectral tests, we next inspected 
the CFH12K images to search for plausible broadband 
counterparts to further confirm that the candidate had 
no strong continuum emission at wavelengths blueward 
of the candidate emission line. This broadband imaging 
test does not absolutely rule out any chance that a given 
candidate is an interloper but instead achieves a more 
limited, though still important goal of verifying that no 
evidence exists that it is one. Finally the three (or more) 
raw DEIMOS spectral frames were inspected at the lo- 
cations of each of the candidates to make sure that the 
candidate emission line was not a cosmic ray artifact. 

3.2. The catalog 



At the end of the search and weeding procedure de- 
scribed in § 3.1, out of an initial sample of several hun- 
dred serendipitous line emitters, we retained 19 single- 
line objects. Each of these nineteen was then assigned 
a confidence class, with class 3 objects being highest- 
quality, and class 1 objects being lowest; in § 4 we do 
further tests on the objects of the different classes to gain 
further confidence in our highest-quality (class 3 and 2) 
objects. 

The confidence classes are defined as follows (see also 
Table 2 for a summary). Class 3 and 2 objects are our 
highest quality LAE candidates — both classes pass all 
the low- z interloper tests described in § 3.1 and are char- 
acterized by a clearly asymmetric line profile typical of 
Lya. Class 3 objects are spatially well-separated from 
foreground galaxies (the primary DEEP2 spectroscopic 
targets or serendipitous foreground objects). Class 2 ob- 
jects are close in projection to a foreground galaxy and 
so their fluxes could potentially be attenuated by pas- 
sage through the intervening galaxies' ISM and/or be 
boosted by galaxy-galaxy lensing. Thus, both class 3 
and 2 objects are very likely to be high-z LAEs, but 
class 2 objects are excluded from certain parts of our 
further analysis because of potential uncertainty in their 
true luminosities. In contrast to class 2 and 3 objects, 
class 1 objects do not show clear asymmetric line pro- 
files, but otherwise are single-line objects that pass all 
the interloper-rejection tests described in § 3.1; they are 
the least secure candidates that we retain in our catalog, 
noting that they probably represent a mixture of LAEs 
and lower- z interlopers (§4). 

The properties of the 19 objects are presented in Ta- 
ble 1, where we list their names (ID), observed wave- 
lengths, redshifts (assuming they are Lya), positions (as- 
suming location is at the centre of the l"-wide slit), limits 
on fluxes and luminosities (at their Lya redshift), confi- 
dence classes, and the redshifts of the primary DEEP2 
targets. Object names are derived from a combination 
of DEEP2 mask number (the first four digits) and slit 
number (the digits that follow the period); the first digit 
of the name identifies the DEEP2 field in which the ob- 
ject is located. Fluxes and luminosities listed in Table 1 
are lower limits because all objects are affected by slit 
losses (which cannot be known for individual objects 
in view of the degeneracy between position and wave- 
length), while, additionally, some objects were observed 
in non-photometric conditions. The DEEP2 spectra were 
not fluxed and so we derived flux calibrations based on 
known instrumental throughputs, which we tested by 
comparing observed broadband magnitudes with spectral 
continua in a larger sample of DEEP2 primary targets. 
On the basis of these comparisons, we believe the LAE 
flux limits reported in Table 1 to be accurate to ^20%. 

In Figs. 2 and 3 we show the direct images and spec- 
tra of our 19 objects. The direct images show 17.4" x 
17.4"regions and are taken from the BRI CFH12K im- 
ages. These images are centered on the positions of the 
primary DEEP2 spectroscopic targets, while positions of 
our line emitters are inferred from the spectroscopic data 
and are marked with open circles. The spectra in Figs. 2 
and 3 show both 2D spectral images and ID spectral 
extractions of the region around the LAE. 

Figure 4 shows the observed fluxes (i.e., lower limits 
on true fluxes) and redshifts of our LAEs. Observed 
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fluxes range from ^5x10 17 erg cm 2 s 1 down to a few 



xlO~ is erg cm" 2 s~ x . Redshifts range over z(Lya)^4.3 
6, as one would expect from the wavelength sensitivity of 
the survey as shown in Fig. 1. Note that there is a con- 
centration of class 3 and 2 objects at z(Lya)=4.4- 4.8. 
This concentration could be a result of real large-scale 
structure, but it could also reflect the larger sky line- free 
volume at these lower redshifts (Fig. 1), combined with 
the fainter luminosities accessible there as compared with 
higher redshifts. Since the bulk of our search was carried 
out in only one of the four DEEP2 fields, we will defer 
the analysis of clustering until the full DEEP2 dataset 
has been searched. 

4. CONFIDENCE TESTS 

The objects in our catalog of 19 galaxies have already 
passed all the individual tests described in § 3.1. In this 
section we present a number of tests that further build 
confidence in the high-z nature of our sample. 

4.1. Averaged spectra 

Averaged spectra give higher signal-to-noise ratios 
than individual ones and thus can give additional confi- 
dence that, as a group, our candidates are genuine LAEs. 
We constructed averaged spectra of three groups of ob- 
jects: (a) Class 3 and 2 objects — i.e., our best LAE can- 
didates, (b) class 1 objects, and (c) 61 serendipitously- 
discovered DEEP2 single-line objects that have earlier 
been ruled out as LAEs through broad-band photometry 
(see Kirby et al. 2007). These are shown in Figs. 6 and 7. 

To produce averaged spectra of each group, we shifted 
the peak of the line in each individual spectrum to a 
common rest-frame wavelength and then averaged each 
group of spectra with inverse variance weighting. We 
employed two types of averaging: uniform weighting 
and line weighting. In the latter case, each spectrum 
is weighted by the signal-to-noise ratio of its emission 
line. Note that the class 1 object 1112.73 was not used 
in making the averaged spectra because this line emitter 
is too close to the principal DEEP2 target galaxy to allow 
for a reliable spectral extraction that's uncontaminatcd 
by flux from the foreground galaxy. 

Line shapes 

Figure 6 shows the averaged (composite) spectra, fo- 
cusing on regions around the discovery line. Also shown 
arc best-fitting Gaussian profiles. The bottom panel 
(panel (c)) of the figure shows the composite spectrum 
of 61 single-line serendipitous DEEP2 objects that have 
been photometrically ruled out as LAEs (Kirby et al. 
2007). This low- z composite spectrum can be used as 
a reference for comparing with the spectra of our LAE 
candidates. The composite spectrum of our nine best 
LAE candidates — i.e., class 3 and 2 objects — is shown 
in Fig. 6(a). It is not well fit by a Gaussian profile and 
instead exhibits the asymmetric line shape expected of 
Lya: the blue side has a steep drop while the red side 
has an extended tail. This asymmetric shape is expected 
given that all the class 3 and 2 objects that make up this 
composite spectrum individually have asymmetric line 
profiles. Overall, this composite spectrum further sup- 
ports the Lya nature of the line in class 3 and 2 objects. 
We present further analysis of this spectrum in § 6. The 



composite spectrum of our class 1 objects, presented in 
Fig. 6(b), does not show asymmetry. This lack of asym- 
metry argues against the class 1 objects being LAEs and 
while it is possible that some of the class 1 objects are 
LAEs, it is very likely that the class 1 sample is contam- 
inated by low-z interlopers. These results do not change 
with uniform rather than line weighting and we conclude 
that class 3 and 2 objects contain appreciable numbers 
of Lya emitters but that there may be significant fore- 
ground contamination among the class 1 objects. 

Presence or absence of other emission lines 

Our experience shows that the most likely single-line 
Lya interlopers are Ha, H/3, [O II] 3727, and [O III] 5007 
emitters (see Martin & Sawicki 2004; Martin et al. 2008). 
Objects with visible Ha often show [N II] or [S II] lines. 
Objects with visible H/3 often show [O III] and vice versa. 
If the discovery line is one of the likely interlopers then 
we might expect to see these other lines in the averaged 
spectra. The averaged spectrum of the 61 non-LAE ob- 
jects shows both [N II] (6548 A and 6583 A) and the two 
[S II] lines (6716 A and 6731 A) when we assume that 
the single line is Ha (see Fig. 7). With uniform weight- 
ing, [N II] is not visible, and the significance of [S II] 
decreases but that line is still detectable. In contrast, 
no such lines are found in the composite spectra of our 
LAE candidates — they are seen in neither the averaged 
spectrum of class 3 and 2 objects nor in that of the class 
1 objects. We also checked the composite spectra of our 
candidates for the presence of [O II] 3727, [O III] 4959, 
5007, [S II] 6716, 6730, Ha, H/3, and H 7 and none of these 
lines were found. The lack of extraneous emission lines 
supports the claim that at least most of our candidates 
are Lya emitters. 

4.2. Redshift distribution 

We next investigate the contamination of our sample of 
high-z LAE candidates by single emission line non-LAE 
galaxies (e.g., Kirby et al. 2007). We compare the red- 
shift distribution of our LAE candidates against a sample 
of 61 single-emission line objects that are known not to 
be LAEs because they are detected in the B and/or R 
bands and clearly lack a spectral break across the emis- 
sion line. We do the comparison not in terms of a tradi- 
tional redshift histogram, but rather in terms of the dis- 
tributions of "effective Lya redshifts" which are defined 
as z e ff,Lya= (A cm /1216 A) — 1. In other words, z e ff,Lya 1S 
the redshift that the object would have if it were a LAE. 
The foreground interlopers in our LAE sample are pre- 
sumably counterparts of the (known) low-z single emis- 
sion line objects, differing from them only in that they 
have a fainter continuum level such that they go unde- 
tected in the B and R bands. 

The bold histogram in Fig. 8 shows the z e ff,Lya distri- 
bution for our sample of 19 class 3, 2, and 1 objects, while 
the hashed histograms are for class 1 (in one case) and 
for class 3 and 2 (in the other) objects only. The light- 
colored histograms show the distribution of the known 
non-LAE objects, with the hashed histogram represent- 
ing the actual distribution and the solid histogram show- 
ing that distribution scaled down by a factor of six. 

As expected, the sample of non-LAE objects is dis- 
tributed more or less uniformly in z e ff,Lya across the 
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available wavelength window. Since the emission line 
in most of these non-LAE objects is expected to be Ha 
(Kirby ct al. 2007), it is not surprising perhaps that 
there are almost no objects with z c gxya < 4.397, where 
[A rcst (Ha)/1216 A] - 1 = (6563 A/1216 A) - 1 = 4.397. 
The rise in numbers from there to Zoff^ya ~ 5 is likely a 
result of increasing volume probed, while the decrease be- 
yond that is likely because of the higher density of night 
sky emission lines at redder wavelengths (Fig. 1) and the 
fall in instrumental efficiency at the longest wavelengths. 

In contrast to the non-Lya objects, the distribution of 
our 19 LAE candidates is peaked at z e ff,Lya ~ 4.7 with 
an extended low-level tail out to the upper end of the 
available range. The distribution of class 3 and 2 objects 
is strongly clustered around z e ff,Lya' v '4.7 and — with the 
exception of object 3307.61 at z=6.0559 — does not show 
a high-Zofl^Lya tail that would be consistent with the low- 
z population. The distribution of class 1 objects is also 
different from that of the non-Lya objects, although the 
class 1 objects do have a more prominent tail to higher 
redshifts. It is thus quite plausible that some fraction of 
the class 1 objects are foreground interlopers. 

Scaling down the distribution of non-LAE objects by 
a factor of 6 causes it to roughly match the extended 
low-level, high-Zofl^Lya tail of the LAE candidate distri- 
bution beyond z e ff,Lyc* ^ 5. The factor of 6 is empirical 
and arbitrary because we have no a priori knowledge of 
the ratio of the number of non-LAE objects with de- 
tectable B- and/or _R-band continuum to the number of 
non-LAE objects in which the continuum is too faint to 
be detectable in our imaging data. If the above factor- 
of-six scaling is taken at face value, it would seem to 
indicate that most of the LAE candidates with z e ff,Lya 
< 5 — including a number of the class 1 objects at these 
z e s,hya — are true LAEs, while the class 1 objects with 
z e s,hya > 5 are mostly single emission line foreground 
galaxy contaminants. 

We conclude that the z e ff,Lya distribution of our LAE 
sample is different from that of the confirmed low- 2; single 
emission line objects, indicating that the LAE candidates 
are drawn from a different population than the low-z 
objects. Given the additional evidence in the form of 
the shapes of the emission lines, we conclude that our 
class 3 and 2 LAE candidates are indeed high- z galaxies, 
while the class 1 candidates may contain a significant 
fraction of low-z interlopers, although it is impossible 
at this point to tell which individual class 1 objects are 
LAEs and which ones are interlopers. 

4.3. Spectroscopic follow-up 

Finally, we have also re-observed one of our candidates, 
namely the class 3 object 4218.94, using DEIMOS at 
lower dispersion and with an exposure time of 3 hours. 
Figure 9 shows the region around the Lya in these fol- 
lowup data. This spectrum clearly shows the asymmetric 
line profile shape expected of Lya, giving us further con- 
fidence that this object is a z=4.53 LAE. 

In summary, we have high confidence that the bulk of 
our class 3 and class 2 objects are genuine, high-z LAEs 
because both individually and as an ensemble they pass 
the full range of tests that we have subjected them to. 
At the same time there is a significant probability that 
the class 1 candidates are contaminated by foreground 



interlopers: their individual and composite line shapes 
are not asymmetric but resemble those of known low-z 
galaxies, but, on the other hand, their redshift distri- 
bution is different from that of single-line low-z objects, 
arguing that a significant number — perhaps around half 
— of class 1 objects could be LAEs; the absence from our 
composite spectra of emission lines expected of low-z in- 
terlopers ( [N II] , [S II] , etc) does not confirm our class 1 
objects as LAEs, but does instead keep them as a viable 
candidate population. 

We conclude that we have high confidence that the 
bulk of our class 3 and 2 objects are genuine high-z LAEs, 
while there is significant probability that the class 1 ob- 
jects are contaminated by foreground interlopers. In the 
following sections we use these samples of objects — re- 
lying mainly on the much higher quality class 3 and 2 ob- 
jects — to constrain the number density of high-z LAEs 
and to study the physical properties in these objects from 
the Lya emission line profile. 

5. NUMBER DENSITY OF LYMAN ALPHA EMITTERS AT 

Z~4.7 

5.1. The simple approach 

A large volume free of OH sky-lines is accessible to us 
in the redshift interval z=4.4-4.9 (Fig. 1) and, in fact, all 
but one of our nine highest-quality (class 3 and 2) objects 
are found in this redshift range (Fig. 4). High-z LAEs are 
a strongly clustered population (e.g., Ouchi et al. 2003), 
and we see signs of clustering in our own data: Fig. 5 
illustrates that our Lya candidates are highly clustered 
on the sky with large parts of the field completely void of 
them. The concentration of redshifts around z^4.7 could 
thus be a result of large-scale structure. However, it is 
also the case that our survey is by far most sensitive to 
finding LAEs around z^4.7 because of the combination 
of the large differential volume, dearth of strong skylines 
at the corresponding wavelengths (Fig. 1), and fainter 
limiting luminosity than is reached at higher redshifts. 
While clustering of LAEs is certainly important in esti- 
mating their luminosity function, we defer the analysis of 
cosmic variance effects until our sample includes fully all 
four of the DEEP2 fields. For now, we take the observed 
numbers of objects at face value and proceed with the LF 
analysis, noting only that our results are based largely 
on a single field (Field 4 of DEEP2), although one that 
spans a large volume. 

The left axis of Fig. 10 shows the cumulative number 
counts of our objects. The open circles show all the class 
3 objects within z(Lya)=4.4-4.9, while the filled circles 
also include the class 2 and 1 objects. The latter are 
therefore likely to be contaminated by foreground inter- 
lopers that possibly lurk among the class 1 objects as well 
as the potentially lensed and/or ISM-attenuated class 2 
objects. The filled circles thus show the more conser- 
vative estimate and henceforth we focus on this sample. 
Note that the DEEP2 points shown in Fig. 10 are lower 
limits on the true number density of LAEs for the fol- 
lowing two reasons. First, our by-eye search no doubt 
missed some objects, particularly at the faint end of the 
population. Second, our objects' luminosities are strictly 
speaking lower limits because lacking detailed positional 
information we have no way of correcting for flux slit 
losses on an object-by-object basis; moreover, some of 
our objects were observed in non-photometric conditions 
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further making their fluxes lower limits. These two ef- 
fects cause our number counts to be lower limits (the first 
limit is in the number direction, the second is in the lu- 
minosity direction) , with the area above and to the right 
of our points in Fig. 10 being permitted by the data. 

To translate from the cumulative number distribution 
for our survey (left-hand axis of Fig. 10) to a cumulative 
number density distribution (right-hand axis) we need 
to divide by the survey volume. Here we use the vol- 
ume that corresponds to integrating over the comoving 
distance from z=4.4 to z=4.9 multiplied by the area 
subtended by the searched slits as a function of red- 
shift (i.e., the integral, from z=4.4 to z=4.9, under the 
lower curve in Fig. 1). This operation gives a volume 
of 2.89 x 10 4 Mpc 3 . Our choice of Az is motivated by 
the fact that almost all of our class 3 and 2 LAEs lie in 
the 0=4.4-4.9 interval, and by the fact that the bulk of 
OH skyline-free survey volume lies in that rcdshift range 
(Fig. 1). Adopting a larger redshift range would decrease 
the number density. In the extreme that allows the full 
(but skyline-reduced) survey volume from the hard lim- 
its of z=4.2 to z=6.6, the volume would increase by a 
factor of 2.4 from the 2.89 x 10 4 Mpc 3 that we adopted 
to 6.92 x 10 4 Mpc 3 ; this effect would result in formally 
lowering the DEEP2 number density shown in Fig. 10 
by a corresponding factor of 2.4. In practice, however, 
the high-z end of this longer rcdshift range would sample 
only the most luminous and rare objects, and thus the 
cumulative number densities are unlikely to be lower by 
as much as that factor of 2.4 from those calculated using 
the more restricted z=4.4-4.9 interval. 

Our assumption about the survey volume here is in- 
accurate for the following reason. Intrinsically luminous 
objects located outside our slits but sufficiently close to 
them could potentially contribute to the observed counts 
by virtue of flux spilling over into the slit; this effect 
could, in principle, lower our number densities. On the 
other hand, intrinsically faint objects will drop out of our 
sample if they are located sufficiently close to the edge of 
the slit to lose significant amounts of their flux due to slit 
losses; this effect acts in the opposite direction, namely to 
increase the true number density. These effects are simi- 
lar to one that affects LF estimates in narrowband imag- 
ing surveys. There, most narrowband filter transmission 
curves are not perfect top-hats but are better approxi- 
mated by Gaussians, and, consequently, in the absence of 
a precise spectroscopic redshift, an object's true flux can- 
not be determined because its position with respect to 
the filter profile is unknown. This uncertainty propagates 
into a bias in luminosity function estimates (further com- 
pounded by any redshift-space clustering that preferen- 
tially favors one part of the filter transmission curve over 
the rest), just as our uncertainty of the object's precise 
position in the slit propagates into a similar uncertainty 
in our simple determination of the LF. The effect can 
be corrected on a statistical basis in both imaging (Shi- 
masaku et al. 2006) and spectroscopic surveys (Martin & 
Sawicki 2004; Martin et al. 2008), but many narrowband 
imaging surveys do not make this statistical correction, 
assuming instead a top-hat filter transmission curve. In 
the present section we take the similar approach of as- 
suming a top-hat slit transmission function, and we defer 
the more accurate approach until § 5.2. 

We compare our straightforward number counts to 



those from the large narrow-band-selected LAE samples 
of Dawson et al. (2008) and of the Subaru Deep Sur- 
vey (SDS; Ouchi ct al. 2003; Shimasaku et al. 2006; 
Kashikawa et al. 2006), which contain some of the largest 
samples of LAEs published to date. We also show the re- 
cent results from the dedicated Magellan z—5.7 multislit 
spectroscopic search by Martin et al. (2008). For the 
Ouchi et al. (2003) results, we show two curves: the up- 
per curve shows the number counts of narrowband-excess 
objects selected using a simple EW cut similar to that of 
Hu et al. (1998); the lower curve is for their sample after 
it has been reduced by a color-color cut correction that 
eliminates objects without a significant Lyman break, as 
well as a statistical correction for foreground interlopers 
based on Monte Carlo simulations (see Ouchi et al. 2003). 
Our adopted redshift interval z— 4.4-4.9 is most di- 
rectly comparable with the z~4.5 results of Dawson et 
al. (2008) and the z-4.9 results of Ouchi et al. (2003). 
The number counts of objects of all three of our confi- 
dence classes are in good agreement with the results of 
these two narrowband imaging surveys. If we consider 
our class 3 objects only then the DEEP2 number densi- 
ties are a factor of a few below those of the narrowband 
imaging surves. This apparent discrepancy is entirely ex- 
pected, however, because our by-eye search is less than 
100% complete and the flux estimates of our LAEs suf- 
fer from (unknown and unknowable) slit losses that may 
make them appear fainter than they really are. Overall, 
we conclude that our serendipitous spectroscopic tech- 
nique is effective at finding samples of high-z LAEs. 

5.2. Accounting for slit flux losses 

In this section wc attempt to estimate the LAE lumi- 
nosity function by properly correcting for flux slit losses 
in a statistical way. We use the statistical LF estimation 
technique developed for this purpose by Martin & Saw- 
icki (2004). The technique is described in detail in Mar- 
tin & Sawicki (2004) and here we review it only briefly. 

The approach is to make a set of parametric models of 
the underlying LAE population and to simulate the ob- 
servational and instrumental processing of these model 
populations to arrive at a predicted number of objects 
for each set of parameters. A comparison of these pre- 
dictions with the observed number of objects (or, in our 
case, lower limit on the number counts, given the intrinsic 
incompleteness of our by-eye selection), then constrains 
the parameters of the underlying luminosity function. 

We assume an LF described by the Schechter (1976) 
function, 
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with freely- varying number density ((f>* ) and characteris- 
tic luminosity (L*) but with a fixed faint-end slope (we 
adopt two values: a=— 1.2 and a=— 1.6). In principle 
the faint-end slope could also be a free parameter, but 
here we have fixed it given that our DEEP2 survey does 
not probe sufficiently deep to constrain the faint end. 

For a given set of Schechter function parameters, we 
expect our survey to find the following number of objects: 
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We set the limiting luminosity L m i n based on the limit- 
ing flux, which we set to / m , n =3xl0 -18 as this contains 
most - but not all of the objects in our catalog. 
Given the subjective, by-eye nature of our search, we do 
not know the survey completeness, £; here we set it equal 
to unity thereby guaranteeing that the expected number 
of objects, N, is a hard upper limit. The survey volume 
V is luminosity-dependent since it decreases for less lu- 
minous objects as intrinsically faint objects need to be 
closer to the middle of the slit to ensure that they have 
sufficient flux after slit losses to make it above the detec- 
tion threshold. As in Martin & Sawicki (2004), we com- 
pute the luminosity-dependent survey volume from the 
cosmological volume increment per unit redshift, the ac- 
cessible redshift interval (which we set as 2=4.4-4.9), the 
effective slit area of the present survey (21.30 arcmin 2 ), 
and the slit loss function defined for our 1" slits and 0.8" 
Gaussian seeing. 

The contours in Fig. 11 show the expected number of 
z=4.4-4.9 LAEs brighter than 3xl0 _18 ergcm~ 2 s _1 for 
different combinations of Schechter <jy* and L* (with a 
fixed at —1.2 and —1.6). There are five class 3 LAEs 
above that flux limit in our survey, which — given that 
our survey is incomplete — gives a hard lower limit on 
the true number of LAEs. Thus, the shaded region of 
parameter space to the upper right of the "5" curve is 
permitted by our observations. Although our current 
catalog is certainly incomplete, it seems unlikely that we 
miss more than nine out of ten true LAEs, and so we set 
a somewhat arbitrary but illustrative upper limit of <50 
LAEs in our plot. The shaded regions in Fig. 11 shows 
the regions of parameter space permitted by these two 
bounds. Note that there difference between the results 
for two faint-end a are very minor — a consequence of 
he fact that at the depths reached by our study we are 
just starting to probe into the faint end of the LAE LF. 

The Schechter fit results reported by a number of re- 
cent Lya surveys are plotted in the bottom panel of 
Fig 11. The results shown are those of Gronwall et al. 
(2007) at 2=3.1, Malhotra & Rhoads (2004) at 2=5.7 
and 6.5, Shimasaku et al. (2006) at 2=5.7, Kashikawa et 
al. (2006) at 2=6.5, and Dawson et al. (2007) at 2=4.5. 
Although different authors assume (or measure) different 
values of a in their fitting, most of these a are close to 
a= —1.6 and, at any rate, the differences in our derived 
LF parameters are not strongly dependent on a. We 
are thus justified in plotting the results of these diverse 
studies together in one common, a=— 1.6, panel without 
explicitly adjusting for differences in their a values. 

When inspecting Fig. 11 one should note that, in con- 
trast to many narrowband imaging surveys, our spectro- 
scopic survey has no equivalent width detection thresh- 
old and therefore in principle, and barring issues to do 
with our survey completeness, is capable of detecting all 
LAEs down to a limiting line flux. In contrast, some of 
the narrowband surveys do impose an equivalent width 
cut, and as a result would miss some of the objects that 
our data arc sensitive to. 

It is reassuring to note that our DEEP2 results are 
in reasonable agreement with Schechter function fits re- 
ported by other LAE surveys at similar and higher red- 
shifts. Overall, then, the agreement between our results 
and those of narrowband imaging surveys that is illus- 
trated in Fig 11 bodes well for the future of multi-slit 



spectroscopic searches as a technique for studying the 
LAE population. 

6. THE SPECTRAL SHAPE OF THE Lya LINE 

The spectral profile of the Lya line carries information 
about the conditions in and around the galaxy that emits 
it. Here we fit the stacked, composite spectrum of our 
nine class 3 and 2 objects (§ 4.1.0) with simple models 
that attempt to model these conditions. Our compos- 
ite spectrum consists of only nine individual spectra of 
relatively poor S/N, so we do not attempt to fit sophis- 
ticated models such as those of Hansen & Oh (2006) or 
Verhamme, Schaerer, & Maselli (2006), but rather use 
simple Gaussian toy models instead. 

One strong limitation of fitting stacked Lya spectra 
which we wish to point out at the outset is that spectral 
averaging relies on shifting the individual spectra to a 
common rest frame. This rest frame can be determined 
only imperfectly for each individual object because of the 
stochastic nature of the IGM absorption of Lya photons 
that truncates the blue wing of the line. Consequently, 
the stacked spectrum suffers some unknown degree of 
"smearing" which cannot be accounted for in a simple 
analysis such as that presented here. Although this effect 
limits the usefulness of the analysis, we press ahead in 
hope of gaining some, even if imperfect, insights. 

The first of our toy models is a single emission line 
that is meant to replicate emitting gas within a high- 
2 galaxy that's broadened by large-scale motions within 
the object. The second model is a two-component model 
that adds a second Lya line redshifted with respect to 
the main line. The physical motivation for the two- 
component model is that of backscattering of Lya pho- 
tons from an outflowing galactic wind or expanding gas 
shell that's powered by the galaxy's star-forming activ- 
ity (e.g., Dawson et al. 2002; Mas-Hesse et al. 2003; Ahn 
2004; Westra et al. 2005; Hansen & Oh, 2006; Kashikawa 
et al. 2006). 

Both models consist of a Gaussian-profile "main" Lya 
line with a width, height, and central wavelength that are 
free parameters of the fit. In the two-component model, 
the second, back-scattering line is likewise modeled us- 
ing a Gaussian profile with a width, height, and central 
wavelength that are free parameters. The velocity offset 
between the main and backscattering lines is also a free 
parameter of the two-component model. In both models 
we mimic the absorbing effect of the intervening IGM by 
setting to zero all flux blueward of the centre of the main 
Lya line. We then convolve the resulting spectra with 
a Gaussian smoothing kernel to simulate the effects of 
the instrumental point spread function (PSF). Since we 
don't know the sizes and in-slit placement of the objects 
that make our composite spectrum, we cannot determine 
the instrumental PSF from the data by using, e.g., night 
sky-lines profiles. Instead, we also make the instrumen- 
tal PSF FWHM a free parameter of the fit. Once we 
have generated a large set of model spectra that reflect 
a large grid of parameter combinations, we then find the 
best-fitting models by means of x 2 minimization. 

Figure 12 shows a comparison of the best-fitting mod- 
els with the data. The observed spectrum is the same 
composite as that shown in Fig. 6(b) but unlike in Fig. 6, 
here the data have not been smoothed. The dashed line 
shows the best-fitting single-component fit. The under- 
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lying (i.e., not truncated) Lya line profile here has a 
FWHM of f .7 A, which corresponds to a velocity disper- 
sion of ct„ = 178 kms -1 . This model, which simulates the 
Lya line as a Gaussian truncated in the blue by inter- 
vening intergalactic gas, accounts well for the observed 
asymmetric Lya line shapes (e.g., Dawson et al. 2002; 
Hu et al. 2004). The width of the line is caused by large- 
scale motions of gas inside the galaxy or possibly by ex- 
tragalactic outflows, while its truncated blue side is due 
to absorption by intervening neutral gas along the line 
of sight. The velocity dispersion of 178 kms -1 in our 
best-fitting single-component model is comparable to the 
200 kms -1 found in the composite Lya spectra at z~5.7 
by Hu et al. (2004). 

A failing of the single-component model is that it un- 
derpredicts the flux at ^1217 A. This excess flux could 
be the signature of back-scattering of Lya photons off of 
a galactic outflow (e.g., Dawson et al. 2002; Mas-Hesse 
et al. 2003; Ahn 2004; Westra et al. 2005; Hansen & Oh, 
2006; Kashikawa et al. 2006). Our two-component model 
aims to simulate this galactic outflow scenario. The best- 
fitting two-component solution is shown in Fig. 12 with 
the thick solid line. It gives a better match to the data 
at ^1217 A than does the single-component model, while 
providing an equally good fit in other regions of the spec- 
trum. In the best-fitting two component model, the pa- 
rameters of the main Lya line are little different from 
those of the single-component model, including a FWHM 
of 1.40 A, which corresponds to er„=147 kms -1 . The 
velocity offset between the main Lya line and the out- 
flow line corresponds to 420 kms -1 — an offset that is 
similar to those found in other Lya emitters (e.g., Daw- 
son et al. 2002; Westra et al. 2005; Kashikawa et al. 
2006). At the same time, the width of the outflowing 
line (cr„=167 kms -1 ) is lower than the 260 kms -1 of 
Kashikawa et al. (2006) at z~6.5, and in the single lowcr- 
z objects analysed by Dawson et al. (2002; 320 kms -1 ) 
and Westra et al. (2005; ^405 kms -1 ). In low-z galax- 
ies, outflow speeds correlate with galactic escape velocity 
(Martin 2005), so these velocity differences may suggest 
that the Lya emitters in our sample are lower-mass ob- 
jects than those studied by the aforementioned authors. 
However, while it is tempting to equate the observed ve- 
locity offsets with outflow speeds, the relation between 
these two quantities is probably less straightforward (see, 
e.g., the models of Hansen & Oh, 2006); unfortunately, 
at the moment our data are not good enough to allow 
more detailed modeling of the properties of the outflows. 

Fig 12 also shows a possible flux excess to the blue of 
the main Lya line. We note that if this excess is real then 
it could be due to some Lya photons escaping absorption 
by the intervening material akin to the escape of ionizing 
continuum radiation (e.g., Shapley et al. 2006; Iwata et 
al. 2008). If real, the escape fraction of Lya photons 
in the blue wing of the line would be ~15%. It will be 
interesting to investigate this issue with the full DEEP2 
LAE dataset. 

Overall, we conclude that the two-component fit that 
attempts to model a backsattering of Lya photons from 
an outflowing galactic wind gives a good match to the 
observed composite spectrum. We regard this result as 
consistent with the idea of large-scale outflows of mate- 
rial in high-z galaxies. 



7. DISCUSSION AND SUMMARY 

Deep LAE surveys are an important tool for the study 
of formation and evolution of galaxies. Despite their 
many very important successes, narrowband imaging 
surveys for LAEs have their limitations: the numbers 
of spectroscopically-confirmed objects are still modest 
and the lack of detailed spectroscopy limits the ability 
to study in detail such properties as galactic outflows. 
Moreover, the lack of spectroscopic redshifts for all the 
LAEs introduces biases in luminosity function and cor- 
relation function studies. Thus, while we must continue 
to press with narrowband imaging surveys, we must also 
develop other, complementary techniques. Martin et al. 
(2008) compare in detail the merits of spectroscopic and 
narrowband imaging imaging LAE surveys and discuss 
the future prospects for these techniques in the coming 
era of extremely large telescopes; here, we focus more 
narrowly on the near-term prospects of exploiting the 
remaining ~80% of the DEEP2 spectroscopic database. 

In this paper we have demonstrated that a large blank- 
sky spectroscopic survey can find significant numbers of 
LAEs. Blank-sky spectroscopic searches have been at- 
tempted in the past (e.g., Crampton & Lilly 1999; Mar- 
tin & Sawicki 2004; Tran et al. 2004; Martin et al. 2008), 
although the number of LAEs fund by them to date 
is very small. Here we have piggy-backed on a large 
intermediate- z survey that is DEEP2 and in just ~20% 
of these data we have found 9 objects that are bona fide 
LAEs (and have a further 10 of which some may also be 
such) . 

So far, ours is a modest sample, but we have been able 
to use it to show that spectroscopic LAE searches can do 
well at constraining properties of the population. Specif- 
ically, we found that the number density of objects in 
our survey is consistent with those found by narrowband 
imaging searches (§ 5), and — using only our stacked 
discovery spectra with no additional follow-up data, we 
confirmed the interpretation that LAEs are associated 
with large-scale outflows of material (§ 6). Much more 
could be done with a bigger sample and with deep follow- 
up spectroscopy. 

A simple extension of the by-eye search described in 
this paper to the remaining ^80% of the DEEP2 data 
can be expected to yield five times more LAEs in total, 
i.e. 42 class 3 and 2 objects (or 88 if class 1 objects are 
also included). However, further improvements in effi- 
ciency over that achieved here are possible. Specifically, 
in the present by-eye search we are sure to have missed 
many objects since the human eye is not an optimal tool 
when confronted with vast amounts of data. An auto- 
mated search of the available dataset should be much 
more efficient. Our first experiments with implementing 
such an automated search are very encouraging and show 
that we can recover the objects found by eye as well as 
find additional, missed line emitters. We speculate that 
with an automated search we should be able to double 
the number of LAEs we find per unit area — as is also 
suggested by a comparison of our number counts with 
those of narrowband imaging surveys (Fig. 10). An au- 
tomated search of DEEP2 would thus give us a sample 
of ^90 class 3 and 2 objects, all with spectroscopic red- 
shifts and with high-resolution spectra suitable for the 
study of the details of gas kinematics. Finally, the spec- 
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tral coverage of DEEP2 gives us the potential to probe 
LAEs all the way up to z^6.6 (Fig. 1), thus allowing us 
to study the evolution of the population from near the 
epoch of reionization to z^4.5 using a single, uniform 
dataset. Another clear advantage of an automated search 
is that it would allow us to accurately calibrate the detec- 
tion efficiency through simulations and thus correct for 
it in our LF analysis. Moreover, analyzing the full four- 
field DEEP2 dataset would also allow us to constrain 
the importance of field-to-field variance. Clustering ef- 
fects appear important even in large fields such as those 
from Subaru/Suprimecam and constraining them using 
the four DEEP2 fields should be a worthwhile endeavor. 
In summary, we have shown that the systematic ex- 
amination of a large spectroscopic database can yield 
significant numbers of faint high-z line emitters which 
are immediately suitable for further studies such as the 
determination of their number density (§ 5) or kinemat- 
ics (§ 6). Automating such a search and applying it to 
an even larger spectroscopic database should be possible 
and should yield a large sample of high-quality high-z 
LAEs suitable for a variety of studies. 
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Table 1. Catalog of Lyo candidates. 



ID 


A(A) a 


^Lya 


R.A. (J2000) 


Decl. (J2000) 


flux b 




T T c 

^ Lyoc 


class d 


7 ■ ° 
^primary 


1112.73 


7293 


4.998 


14:16:18.23 


+52:17:05.5 


8.1x10" 


-18 


2.1xl0 42 


1 


0.7464 


1143.85 


7505 


5.1835 


14:15:20.55 


+52:07:38.6 


7.8x10" 


-19 


2.2xl0 41 


1 


0.8789 


1150.116 


7019 


4.7694 


14:16:07.32 


+52:25:40.4 


5.4x10" 


-18 


1.3xl0 42 


3 


0.7777 


2103.71 


7029 


4.7787 


16:47:01.01 


+34:49:56.7 


7.8x10" 


-18 


1.8xl0 42 


3 


0.9727 


3307.61 


8603 


6.0559 


23:32:57.05 


+00:01:24.3 


9.4x10" 


-18 


3.8 xlO 42 


2 


0.8718 


4107.130 


6972 


4.7281 


02:27:28.03 


+00:35:34.3 


9.1x10" 


-18 


2.1xl0 42 


3 


0.8311 


4110.123 


6563 


4.3975 


02:27:59.03 


+00:34:31.1 


1.5x10" 


-17 


3.0xl0 42 


2 


0.9101 


4112.49 


6668 


4.4842 


02:27:58.09 


+00:26:17.3 


1.2x10" 


-17 


2.4xl0 42 


2 


0.9071 


4147.18 


6777 


4.5750 


02:27:25.75 


+00:36:50.7 


6.6x10" 


-18 


1.4xl0 42 


2 


0.9352 


4180.132 


7682 


5.3166 


02:26:08.48 


+00:42:32.0 


5.5x10- 


-18 


1.7xl0 42 


1 


1.0393 


4218.94 


6724 


4.5252 


02:30:48.65 


+00:31:38.2 


1.5x10" 


-17 


3.2 xlO 42 


3 


0.8473 


4240.64 


6769 


4.5670 


02:28:38.04 


+00:41:07.4 


4.9x10- 


-18 


l.OxlO 42 


1 


0.9691 


4243.121 


6499 


4.3452 


02:28:59.67 


+00:47:30.7 


6.5x10- 


-18 


1.2xl0 42 


1 


0.9947 


4243.126 


6996 


4.7518 


02:29:01.06 


+00:48:17.2 


7.2x10- 


-18 


1.7xl0 42 


1 


1.0213 


4256.30 


6919 


4.6883 


02:31:14.91 


+00:37:31.2 


3.9x10- 


-18 


8.8xl0 41 


1 


— 


4257.9 


7415 


5.1026 


02:30:59.78 


+00:47:56.3 


2.6x10- 


-18 


7.0xl0 41 


1 


0.9263 


4260.25 


7908 


5.5019 


02:31:14.92 


+00:37:36.1 


4.9x10- 


-18 


1.6xl0 42 


1 


1.0370 


4280.73 


6854 


4.6359 


02:28:42.80 


+00:35:32.8 


3.3x10- 


-18 


7.3xl0 41 


3 


0.9578 


4280.76 


6865 


4.6483 


02:28:40.70 


+00:35:56.8 


1.1x10" 


-17 


2.4xl0 42 


1 


1.0045 



a Wavclcngth of the peak of the emission line 

b In units of erg/s/cm 2 and not accounting for the (unknown) slit losses 
c In units of erg/s and not accounting for the (unknown) slit losses 

d Quality: 3=best: asymmetric line shape; 2=asymmetric line shape but object is very close to a foreground galaxy; 1= 
line consistent with Lya, but without clear asymmetric line shape. 
c Redshift of the primary spectroscopic target 



isolated emission 



Table 2. LAE candidate confidence classes 



class 



description 



number of objects 



3 isolated object with asymmetric line profile characteristic of Lya; passes all the tests in 

2 as class 1, but is close to foreground galaxy and so its flux could be inaccurate 

1 line profile not clearly asymmetric, but otherwise object passes all the tests in § 3.1 



3.1; our best candidates 
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Fig. 1. — The wavelength coverage of DEEP2, and the survey's volume as a LAE survey. The top panel shows the volume surveyed 
in the present work which is based on 83 masks (lower thick line) and that covered by the full DEEP2 survey of 385 masks (upper thick 
line). The volume shown is that subtended by the full DEEP2 slits — i.e., no correction for obscuration by low-z galaxies or slit-end effects 
have been applied. Note that the volume is shown per unit redshift (Az=l). The vertical gray lines mark regions of high sky noise due 
to atmospheric emission lines. The lower panel shows the cumulative volume of the present survey, with the dashed curve showing the 
volume without correcting for wavelength regions made unusable by night skylines and the solid curve showing the volume after applying 
such a correction. While there are many skylines in the red, DEEP2's spectral resolution provides enough clear inter-skyline regions to give 
significant volume for LAE searches at higher redshifts. 



APPENDIX 
NOTES ON INDIVIDUAL OBJECTS 

All of the candidates discussed below have no associated spectral continuum either on the red or blue side unless 
otherwise specified. 

Object 1112.73. This is a poor candidate due to the fact that the candidate's emission line is not asymmetric; its 
location in the slit is very close to the primary target galaxy. The candidate, however, has no other associated lines. 
The object may have a faint counterpart in the B,R, & I images which could rule out Lya as line ID. We retain the 
object as a class 1 candidate. 

Object 1143.85 has a very weak emission line at 7505 A which shows little to no asymmetry in the 2D and ID 
spectra. The candidate is fairly well spatially separated from the primary target and is absent in the BRI images. 
The faintness of the emission line gives the candidate less credibility as a LAE as it could possibly be [O III] at 5007 A 
with undetected 4959 A and 4861 A. We classify the object as a class 1 candidate. 

Object 1150.116 has a strong single emission line at 7020A that is moderately asymmetric in both the ID and 2D 
spectra. The object is located well away from the primary target galaxy. There may be some continuum associated 
with the candidate on the redward side, especially a few hundred Angstroms redward of the potential Lya line. The 
candidate is undetected in the BRI imaging. We classify it as a class 3 object. 

Object 2103.71 has a very strong emission line at 7030 A that is strongly asymmetric in both the ID and 2D spectra. 
The object is located far from the primary DEEP2 target and is undetected in the BRI images. We classify it as a 
class 3 object. 

Object 3307.61 has an asymmetric emission line at 8603 A in both the ID and 2D spectra. Unfortunately, the 
object is very close to the primary target making it difficult to extract without contamination from the primary 
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Fig. 2. — Direct images and spectra of our objects of confidence class 3 and 2. The three left-most panels show direct B, R, and I 
CFH12K images centered on the primary DEEP2 spectroscopic targets. Slit positions are indicated as are the inferred positions of our 
Lya candidates. These broadband images are 17.4" X 17.4" CFH12K image segments. The right-most panels show the spectral images of 
the wavelength regions around the Lya candidates; wavelength increases to the right and the length of the image corresponds to 33 A for 
each object, while the height of the image shows the full extent of the spectrum along the slit. The 1-D spectra show the line shape at the 
position of the candidate line, with the wavelength scale in angstroms. 



target continuum and thus we cannot search for spectral continuum of this object. The object is undetected in the 
BRI images, but the location of the object is so close to the primary DEEP2 target that even a bright broadband 
counterpart would be difficult to detect. We classify this object as a class 2 object. 

Object 4107.130 has a strong emission line at 6973 A that is moderately asymmetric in both the ID and 2D spectra. 
The candidate is not visible in B or R, but could be present in the I image at very low level of significance. Class 3 
object. 

Object 4110.123 has a very strong emission line located at 6564 A that is well separated from the primary target 
galaxy but is nearly superimposed on another serendipitous continuum source. This other serendipitous object has an 
[O II] emission line at 6843 A and so is at z— 0.836. Our candidate 6564 A line does not correspond to any known lines 
at z=0.836. Moreover, the candidate Lya line shows strong asymmetry in both the ID and 2D spectra. The candidate 
appears absent in the B image and is present in both R and I bands, although this is difficult to tell with certainty 
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Fig. 3. — As for Fig. 2 but for the less secure objects of confidence class 1. 



due to its proximity to the serendipitous [O II] emitter. Flux measurement of this LAE candidate could potentially be 
affected by light from the [O II] emitter and hence it is listed as class 2. 

Object 4112.49 has a very strong emission line located at 6669 A that is nearly superimposed on the primary DEEP2 
target galaxy spectrum. The primary target is at a redshift of 0.8091 as determined from very strong [O II] emission 
at 6741 A. The LAE candidate line is not matched by any known features at the redshift of the primary target. The 
candidate Lya line shows strong asymmetry in both the ID and 2D spectra. Broadband images are not useful in this 
case due to the proximity to the primary DEEP2 target. This is a class 2 candidate due to its proximity to the primary 
target. 

Object 4147-18 has a strong emission line at 6778 A that is strongly asymmetric in both the ID and 2D spectra. The 
candidate is also present in both the R and I images but absent in B-band, lending further credence to this candidate's 
bid as a LAE. The potential Lya line is nearly superimposed on another serendipitous spectrum which is an [O II] 
emitter at a redshift of 1.21; however the candidate LAE line does not match any lines at this redshift and hence we 
conclude that the object is not associated with the z—1.21 emitter. While this is one of our best candidates based on 
its spectral shape, its apparent brightness may be affected by the foreground [O II]-emitting galaxies and hence it is 
classified as a class 2 object. 

Object 4180.132 has a very strong emission line located at 7683 A that shows no asymmetry in either the ID or 2D 
spectra, looking very much like a typical H/3 or [O III] line. The position of the emission line is spatially well separated 



C\2 



s 



CD 



X 
3 



2x10 1? 



Lyman a Emitters in DEEP2 

z(Lya) 
4 4.5 5 5.5 6 



15 



10 



-17 




6000 



7000 



8000 



1 , 1 , , 

1 


III 

1 


III 

1 


1 1 1 1 1 1 1 


i 




• class 3 


1 






o class 2 




- 

o 


• 




o class 1 




1 
1 
1 
1 


o 
o 




- 


1 


• 


o 




1 

1 

o 


• 
o 
o 


o 


" 


1 


o • 


o 
o 




1 


o 




- 


1 

ill 1 


. 1 


o 

1 1 1 1 1 1 1 1 


i 



9000 



wavelength(A) 



Fig. 4. — Observed fluxes and wavelengths of our line emitters. The fluxes shown are lower limits on the true fluxes as they have not 
been corrected for slit losses. The horizontal dashed line denotes our adopted limiting flux, and the vertical lines bracket the redshift range 
used in our number density calculations (§ 5). Filled symbols show objects with the two highest Lye? confidence classes. 
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Fig. 5. — Positions and Lya redshifts of our LAE candidates in the 4th DEEP2 field. Small points show the positions of the 11,250 
DEEP2 spectroscopic targets in this field and illustrate the topology of the survey. It is clear that our LAE candidates are clustered both 
is space and in redshift. 
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Fig. 6. — The Lya candidate line shape for three averaged groups of spectra: (a) class 3 and 2 objects, (b) class 1 objects, (c) the 61 
single-line DEEP2 objects that are known not to be LAEs. Dotted lines show Gaussian fits. 



from the primary target galaxy but does not appear to be present in any of the broadband images. Object classified 
as class 1. 

Object 4218.94 has a very strong, very asymmetric line at 6724A that is well separated from the primary target. 
The emission line has a distinct asymmetric shape that is truncated on the blue end and extends in the red. The 
broadband images further lend credence to the supposition that this is a LAE, as the candidate appears in both the 
R and I bands but is absent in the B band. This is quite possibly our best candidate. Class 3. 

Object 4240.64 contains a moderately bright emission line at 6769 A that exhibits little to no asymmetry. This 
emission is fairly close to the primary DEEP2 target spatially, such that some of the target's continuum could have 
corrupted the spectral extraction of the LAE candidate. The candidate appears to be absent in all the broad band 
images. Class 1 object. 

Object 4243.121 has a very bright emission line at 6500 A that shows little to no asymmetry. The broadband images 
show a galaxy at the position of this candidate that is present in R and / but absent in B. However, it is difficult to 
tell whether this galaxy fell within the slit and thus whether it is associated with the line emitter in question. Class 1 
object. 

Object 4243.126 has a moderately bright emission line at 6997 A. The analysis of this candidate is difficult because a 
line of bad pixels runs directly through the emission line's center. However, it appears that the candidate's line shape 
has no asymmetry either in the ID or 2D spectrum. The candidate is undetected in the B band and may be present 
in the R and / bands. Class 1 object. 

Object 4256.30 has a bright emission line at 6919 A that shows no asymmetry and looks similar to a normal Ha or 
[O III] emission line. The candidate is quite close to the primary DEEP2 target spatially, making it challenging to 
tell whether the candidate is visible in the BRI images. However, it does appear absent from these images. Class 1 
object. 

Object 4257.9 has a weak emission line at 7415 A which has no asymmetry. The object does not appear to be present 
in any of the 3 broadband images. However, there may be an object in R & I near the location of the candidate LAE 
and if the LAE candidate is associated with this galaxy then it would be ruled out as a LAE. Class 1 object. 

Object 4260.25 has a moderately bright emission line located at 7909 A which has no asymmetry but — given a 
non-detected continuum — has a large equivalent width that is uncharacteristic of low- 2: emission lines. The ID 
spectrum is affected by bad night sky subtraction making it difficult to tell if there is any continuum associated with 
the emission line. The candidate is not present in cither the B or R broadband images. There may be a galaxy present 
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Fig. 7. — Spectral details for the co-added groups of spectra. The top panels are for the composite of the 9 best candidates (class 3 and 
2 objects), while the bottom panels are for the 61 non-LAEs single-line objects. We assumed that we have detected Ha and not Ly« and 
then use dashed lines to show the expected positions of [N II] (left panels) and [S II] (right panels). 



at the location of the candidate in the I image, but this is not certain. Class 1 object. 

Object 4280. 13 has a strong emission line at 6855A that shows very well defined asymmetry in both the ID and 2D 
spectra. The broadband images show no counterpart in B, R, or /. Class 3 object. 

Object 4280. 16 has a strong emission line at 6866A which has little to no asymmetry but does have a large equivalent 
width making it unlikely that it is a line associated with a lower-z galaxy. The candidate is absent from all three 
broadband images. Class 1 object. 
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Fig. 11. — Constraints on the luminosity function parameters. The expected number of objects brighter than 3x10 erg s _1 as a 
function of the Schcchtcr parameters of the LF. The two panels assume different fixed faint-end a, but, as can be seen from the small 
differences between the two panels, the choice of a is not dramatically important. The shaded region shows the values of the Schechter 
parameters that are permitted by the number of objects we detected (the upper limit of this region assumes that we are missing no more 
than 9 out of 10 objects). In the bottom panel we show the Schechter fit parameters reported by a number of recent Lya surveys: Gronwall 
ct al. (2007, G07), Malhotra & Rhoads (2004, MR04), Shimasaku et al. (2006, S06), Kashikawa et al. (2006, K06), and Dawson et al. (2007, 
D07). 
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Fig. 12. — Model fits of the composite spectrum of our class 3 and 2 objects. The thin solid lines show the two components of the 
two-component model as they would appear in the absence of intcrgalactic absorption. 



